EPJ manuscript No. 

(will be inserted by the editor) 



The effect of clusterings on the equilibrium states of local 
majority-rule: Occurrence probability and robustness 

<N ■ 

Yu-Pin Luo 1 Jinn- Wen Wu 2 and Ming-Chang Huang 3a 
, 1 Department of Electronic Engineering, National Formosa University, Huwei, 63201 Taiwan 



o 



in 



c 

o 
o 



^v^j ' 2 Department of Applied Mathematics, Chung- Yuan Christian University, Chungli, 32023 Taiwan 
3 Department of Physics, Chung- Yuan Christian University, Chungli, 32023 Taiwan 



Received: date / Revised version: date 
Abstract. Abstract is missing. 

PACS. XX.XX.XX No PACS code given 

1 Introduction ity, here we use artificially constructed networks to define 

the neighbors of an unit as its nearestly connected nodes. 



t-H ■ The collective behaviors, such as the emergence of system- 

> : 

, wide coordination in nature, the appearance of the con- 



Because of the locality in LMR, one may expect that the 
distribution of cliques of a system can affect the occur- 



sensuses of opinions in social systems, and other related 



rence probabdity of collective behavior, and this work is 

' phenomena, arc of great interested for many researchers [II , , , , 

■ devoted to analyze such effect. 

^ '[![3lH[5]. The local majority-rule (LMR) has often been 

^ > . employed to study the arising of such behaviors. The rule The global topology of a network can be characterized 

^ ! is simple, based on the principle of majority vote with- by two quantities, the degree distribution and the clus- 

out much consideration in psychological level; it dictates tering coefficient 0. The total number of connections of a 

that the time evolution of the state of an unit (individual node is referred as the degree of the node k, and the prob- 

or agent) is determined by the majority-favored state of ability that a randomly chosen node has k connections is 

its neighbors [5] . The neighbors of an unit can be given g i ve n by the degree distribution P (k); the tightness of a 

by geographic, cultural, social, or organizational proxim- c ii que formed by a site and its directly connected neigh- 



email: mchuang@cycu.edu.tw bors can be characterized globally by the clustering coef- 
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ficient of a network C . The question concerning with the of cliques in the system affects the corresponding trajec- 
role of network topology on dynamical cooperative behav- tory strongly, and hence affects the type of equilibrium 
ior were discussed by Sood and RcndcrJSj and by Suchecki states reached by the trajectory. One of the questions we 
et al.[9] in the voter model which may be viewed as a sta- intend to address in this work is as follows: Starting with 
tistical model of LMR[W\. The mean time for reaching strongly disorder states, what role does the clustering co- 
the state of collective behavior was shown to have differ- efficient of the system play in the type of the equilibrium 
ent scaling behaviors with respect to the number of nodes states reached by the system? 
for different decay exponents 7 of scale-free (SF) networks 

Another question that attracts our attention is the re- 

betwecn 7 > 3, 7 = 3, 2 < 7 < 3, 7 = 2, and 7 < 20. 

lation between the robustness of a equilibrium state and 

Moreover, network geometry was also shown to have im- 

the clustering coefficient of a system. The study in this 

portant effect on the dynamics such as the average survival 

aspect may not only reveal the stability of an equilib- 

timc of metastable states in finite networks, the linear size 

rium state but also provide an estimation for the external 

scaling law of the survival time, and the size of an ordered 

strength required to break the state. The robustness of 

domain[9]. As LMR is the root of the voter model, these 

an equilibrium state can be characterized by its escape 

features may be traced to the properties of the equilibrium 

rate after introducing fluctuation to perturb the system. 

states of LMR. Hence the study within LMR may provide 

In fact, fluctuation is an unavoidable component for real 

more insights to the question. 

systems. The attempt was made by Moreira et al. to in- 
Different 7 values of SF networks characterize the dif- elude noise into the dynamics by changing the transition 
ference in the appearance of hub-nodes. Here the hub- probability of LMR from 1 to 1 — 77, where the param- 
nodes are referred to those possessing large degree of con- etcr r\ characterizes the average effect of fluctuation [T3]. 
nections. The existence of hub-nodes may affect strongly The authors showed that the presence of fluctuation may 
the efficiency of reaching an equilibrium state of LMR. increase the probability and the efficiency of occurring col- 
For this aspect, Zhou and Lipowsky showed that there lective behavior for systems with small-world characters, 
exists categorical difference between the SF networks of In this work, we take a microscopic approach by propos- 
7 < 5/2 and those of 7 > 5/2 for the scaling behavior ing a stochastic LMR in which, each node-state contains 
of the relaxation time from a strongly disorder state to- a component of white noise. As the equilibrium states of 
wards an order state|13j. But there are other equilibrium LMR become transient in stochastic dynamics, this pro- 
states associated with LMR that are different character- posal yields the Arrhenius equation for the escape rates, 
istically from the state of collective behavior. Moreover, Then, we determine the dependence of the prefactor and 
as a system evolves from an initial state, the distribution the activation energy of the Arrhenius equation with re- 
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spect to the clustering coefficient. This information allows A. Here the matrix A is symmetric with the elements a%j = 

us to show the effect of the clustering coefficient on the 1 for the connected sites i and j, and otherwise. The 

robustness of a state explicitly. As the state of collective dynamic variable associated with a site i is denoted as 

behavior is of great interested, we also study the mean X{, which takes two possible values, either 1 or — 1. The 

first-passage time from a strongly disorder state to the system evolves from an initial to a new configuration in 

state of collective behavior. Such study, in addition to the discrete time step according to LMR whose operation can 

robustness of the state, may provide further understand- be either synchronous or asynchronous. In this work, we 

ings about the role of fluctuation in the process of reaching consider the synchronous dynamics for which, the rule can 

the state of collective behavior. be written as 

This paper is organized as follows. In Sec. II we de- / l 

x t (t + 1) = sgn 2^ aijXj (t) (1) 

fine the LMR, classify its equilibrium states, and briefly \i=i / 



describe the generating processes for the networks used in 
the numerical study. Based on the LMR, we numerically 
calculate the occurrence probabilities of different classes of 
equilibrium states for the systems starting with strongly 
disorder states, and the results as functions of clustering 
coefficient are shown in Sec. Ill where the dependence on 
the system size for the occurrence probability is also dis- 
cussed. In Sec. IV, we first introduce the stochastic LMR, ^ = _ x . (fj 



for i = 1,...,N, where the sgn function is a standard 
threshold function with sgn(:r) = +1 for x > and —1 for 



x < 0, and we set Xj (t + 1) = Xi (t) for ciijXj (t) = 

0. The dynamics of Eq. (TTJ) has been widely studied in 
discrete neural networks, and the existence of equilibrium 
states can be shown by employing the Lyapunov energy 
function [5UT5]. 



N 



then the analysis on the escape rates of different classes 
of equilibrium states based on the stochastic dynamics 
are given. Moreover, the results for the mean first-passage 
time to the state of system-wide coordination are also pre- 
sented in this section. Finally, a summary of the results 
and some general conclusions are given in Sec. V. 



N 

X! (* ~ x ) 



(2) 



Moreover, the period of an equilibrium state is either 1 or 

2B51- 

The number of equilibrium states indicates the capac- 
ity of a neural network. But, we are interested in the oc- 
currence of the state of collective behavior and the global 
characters of other equilibrium states in case that the col- 
lective behavior can not be reached by systems. Thus, the 
equilibrium states of period- 1 arc divided into two classes, 
We first specify the LMR and classify the corresponding Sq and Si . Here the class So is for the states of collective 



2 Deterministic Dynamics and Networks 



equilibrium states. Consider a network system with the behavior for which, all node-states have the same value, 
distribution of edges given by an N x N adjacency matrix either 1 or — 1; and the class S\ consists of all trapped 
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states. For the equilibrium states of period-2, the system C values are employed for numerical study in this work, 
oscillates between a pair of states, and we refer all pairs of and they were generated by using the modified schemes of 
equilibrium states as the class 5*2. As the order parameter preferential attachment, proposed by Holme et al.|17j and 
is defined as M s = ^iLi x %J^i the two states of So have by Leary et al.[18], in a systematic way. For tuning the 
the value 1 and —1, respectively, and the states of S\ and C value without altering the 7 value, a step, called triad- 
S2 have \M 8 \ < 1. formation, is added to the process of preferential attach- 
Two types of networks, Watts- Strogatz ( WS) and SF ment with a assigned probability; the larger the assigned 
networks, are used to define the adjacency matrix A in the probability of performing triad-formation, the larger the 
dynamics of Eq. JT]). The WS networks are made from C value of the resultant network is[T7]. Alternatively, we 
a regular lattice for which, N sites are placed around a alter the 7 value without affecting the C value by switch- 
circle and each site has degree, say ko, connecting to the ing the uniform distribution for the random numbers used 
right and to the left symmetrically, then a probability p is in the process of preferential attachment to the distribu- 
assigned to rewire the edges randomly [TBJ. As the p value tion of a designed probability density function [IS]. As the 
increases from to 1, the resultant network changes from designed probability density function further enhances the 
a regular lattice to a random graph with the clustering probability of connecting a new edge to a site with large 
coefficient C decreasing from the highest value down to degree, we obtain a network with a larger 7 value. On 
k /N . Here the C value of a network is defined as the the other hand, for the opposite tendency in the designed 
average of the clustering coefficients associated with all function the resultant network has a smaller 7 value, 
sites, and the clustering coefficient of a site, say i, is given 
as 

3 Occurrence Probabilities of Equilibrium 

d = 2yt (3) 

Mfci_1) ' States 

where is the number of existent edges between the hi 
neighbors of the site i. 

The SF networks are a special category of networks We first calculate numerically the occurrence probabil- 

for which, the degree distributions take the form of power ity for three different classes of equilibrium states of Eq. 

low as P (k) ~ k~ 1 with decay exponent 7. A conventional ([!]) to analyze the clustering effect. The occurrence prob- 

way of generating a SF network is the scheme of prefer- abilities denoted as Pj for the class S% of the equilibrium 

cntial attachment proposed by Barabasi and Albert [TT]. states with i = 0, 1, or 2 for systems modelled by the 

and it yields a network with 7 sa 2.9 and small C value WS and the SF networks. To minimize statistical errors 

C w jy-° ' 5 [X2]. The SF networks with different 7 and in the simulation results, we generate 1000 samples for the 
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WS networks with a given value of rewiring probability, D , AT . no 0.95 + (1.04 x 10 5 /7V 2 07 ) z (N,C) 

P2 (Jy, C) — 1.03 



and the corresponding clustering coefficient C is defined 
as the average of the coefficients of all samples. Then, the 
Pi value associated with a C is given by the fractional per- 
centage of occurrence for the equilibrium states belonging 
to the class Si over the 10 6 trajectories. Here the trajecto- 
ries all start from strongly disorder states and are equally z (N, C) 
distributed in the 1000 samples. 

The results of P; for the WS networks arc shown as 



1 + z (N, C) 

(5) 

and 

Pi (N, C) = 1 - P (N, C) - P 2 (N, C) , (6) 
where z (iV, C) is 

' V. I 5 . (7) 

(0.65/iV - O7 -0.25)J v ; 

As shown as the solid lines in Fig. 1(b), the expressions 
agree with the numerical data very well. Eq. Q reveals an 



the plots of Pi vs. C in Fig. 1(a) for different fc values important fcaturc about p 0; that iB> thc Pq valuc is sup _ 

with N = 1000 and in Fig. 1(b) for different N values pregsed &g increasing N or q. Consequently, there exists a 

with fc = 8. Here the C values arc in thc range < C < gite number N * (fco) that Po ~ for iV > N* (fc ), where 

0.5 for which, the corresponding rewiring probabilities are N * ^ dccrcases as ko deC reases. Moreover, Eq. © indi- 

ranged between 0.04 < p < 0.5; this is also the range of categ that the ^ yalue increages as N increases . ThuSj we 

C generated for the SF networks. As shown in the plots, m&y haye Pq w q> ^ M Q ^ &nd p 2 ^ x ioY N > N * (fco) . 

systems are more likely to be led to the states of collective Thig ig demonstratcd by the numeric al results shown in 

behavior for small C values (C ~ 0.05) , to the equilibrium pig 2 where the p% &g functions of c for TV = 6000 

states of Si for medium C values (C ~ 0.15), and to the 1 , 

1 V / ' and fco = 4 are given. 

equilibrium states of S 2 for high C values (C ~ 0.3) . The Thc statistics for thc numcrical study m the SF net- 
results also indicate that the P value is suppressed and workg ig gimilar tQ the c&ge of ^ m networks _ As the 
the P 2 value is enhanced as the node number N increases, modificd schemcs of prcfcrcntia i attachment are applied, 

and thc tendency is opposite for increasing the average , , ,• , r~, , rr^i 

J * 1 ° b thc 7 value may change slightly m tuning the G value [17], 

number of degree k . and the c yalue may alter slightly m tuning t h e 7 value [E] . 

Analytic expressions Pi (JV, C) , which can fit the data Thus, the samples are characterized by the set- values (7, C) 

shown in Fig. 1(b) properly, are found for a better un- with 1000 members belonging to a (7, C) . There are three 

derstanding about the TV-dependence of the occurrence distinct 7 values, 7 = 2.82, 2.55, and 2.29, and the possible 

probabilities P l . The results are C values for a 7 value are in the range between and 0.5. 

All samples have the site number N = 5000 and the aver- 



1.11 - (0.06) TV - 16 

f°- 35 ]C-6.70}' a&CU ' 

(4) thc 10 6 trajectories for which, they are equally distributed 



Po (N, C) — ; = r- , age degree of a site (k) = 4. The P, value is the result over 

V 7 1 + exp {[91.53- (407.28) //V°- 35 ]C- 6.70} & B W 
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in the 1000 samples of the SF networks belonging to the where (t) is the Gaussian white noise with the zero 

same (7, C) and each trajectory starts from a strongly dis- mean and the ^-function correlation, i,e, (t)) = and 

order state. The numerical results are shown as the plots ^ M£j ( s )) = — s )j an d D is the diffusion con- 

of Pi vs. C in Fig. 3 for three different 7 values. stant which characterizes the strength of noise. The equi- 

The results obtained from the SF networks indicate Librium states of Eq. (JTJ become transient for the stochas- 

that although the 7 value may affect significantly the con- tic dynamics of Eq. © . Consequently, the escape time for 

vergent speed of the system leading to a equilibrium state, different classes of the equilibrium states can be measured, 

it has little effect on the Pi value for which, the C value the results are denoted as (t^) with i = 0, 1, and 2 for 

plays a major role. Moreover, similar to the case of the the class So, Si, and S2, respectively, and the inverse of 

WS networks, as shown in Fig. 3 the clustering coefficient ( T ^) yields the escape rate (m), where the bracket of r"' 

C drives the system from the state of collective behavior or Ki represents the average of the results over the samples 

at low C (C < O.l) to the the phase of oscillation between of different networks. 

two states at high C (0.5 > C > 0.3), and the system is As the stochastic dynamics of Eq. © is applied to 

trapped in a state of Si at medium C. One may further the equilibrium states of Eq. (TTJ) , one can expect that the 

expect that the iV-dependence for Pj has the same feature escape rates obey the Arrhenius equation j!9j. 

qualitatively as that for the WS networks. / \ _ a ( 1 fa\ 

\Ki) — A-i CXp I — I , (J) 

where Ai is the prefactor, and AGi is the activation en- 

4 Stochastic Dynamics ergy Note tnat based on the fluctuation-dissipation the- 



orem we may identify D = ksT with the Boltzman con- 
stant ks and the absolute temperature T[5D]. The factor 
Ai, which is equivalent to the rate constant in a chemical 
reaction, signifies the entropy effect, and one can expect 
that it has a strong dependence on the clustering coeffi- 
cient. The value of AGi gives the maximum potential en- 
ergy required to escape from the equilibrium state. Thus, 
the results of Ai and AGi provide insights on the robust- 
ness of the equilibrium states of different classes as the 
geometric structure of a system varies. 
(t + 1) = s g n ( ~^^o,ij s S n (xj (t) + \Z2D£j (t)J J , We perform the numerical measurements for (Ki) in 

(8) the SF networks, and the statistics for the measurements 



Noise is a very natural component for real systems, its 
physical origins can be traced to incomplete information, 
processing errors, or other environmental perturbations. 
To add a component of noises to the LMR, we propose 
a stochastic version based on the assumption that the ef- 
fect of noise is localized and appears as the fluctuation 
in recognizing the value of a node-state by its connected 
neighbors. Then, the stochastic LMR is given as 

N 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 7 

is described as follows. We first single out the mo equilib- value, and their C-dependcncc is shown explicitly in Fig. 

rium states belonging to the class S, for a set-value (7, G) 5 for 7 = 2.82. 

with mn = 25000, where the mn states of the class So are . , c c am* 1 j r 

' Some important features for Fq. ([9| are revealed from 

identical as a%- = +1 for all nodes. Then, the number of ,, n t^. . , _ t^. ,, , 

' the results shown in .bigs. 4 and 5. firstly, the decay rate 

time-steps required to escape from each of the mn states is r , , , c n v. u , . 

^ n ^ for the state of collective behavior is universal m the sense 

measured with a preassigned cut-off time-step L a „ such a ,, , ., . . , , r , 

r that it is independent of the network geometry, so are 

measurement is repeated for 20 times to realize the Gaus- , r . , , 

r the ZiGi values for the equilibrium states belonging to 

sian distribution of the noise £„■ (i), and the escape time- ,11010™ c 

Sl w ' 1 the classes Si and 62. Moreover, the activation energy for 

step (t^) for the class & is given by the average value , , c a j o • 1 i.i n , r ,1 , 

^ \ / the states of Si and S2 is larger than that for the states 

over the 20 simulations and over the mn states. c o ,l rr ,1 1 , ,1 7=f 

of 00. r or the entropy effect on the decay rate, the G 

Our results for the SF networks with (j,C) = (2.82, 0.302^ lue haS a si S nificant role in determining the values of 

, r . „ , , £1 / a)\ 1 /n -^l an d ^2- As the results of Fig. 5 indicate, the A2 value 

are shown m lig. 3 as the plots of log \t { ') vs. f / D with ' 

n t jr, j,i r j-rr , /— 7=<\ i i decreases for increasing the C value, and the Ai value 
t = 0, f, and 2, and those for different (7, G] values also ° ' 

.1, 1 . . n , mi . ,. behaves oppositely; moreover, the equilibrium states of 

yield straight lines m the same plots. Ihe results indicate ^ J ' 

, , , , , > c . , 1 ft possess a very large prefactor, and this overcomes the 

that the escape rates [Hi) for the equilibrium states be- 

1 • , ,1 i o -.i-v, -C (fm v. a higher activation energy and renders the states to be very 

longing to the class bi agree with Fq. (|9|) , such agreement ° ° J J 

rr,. , r -irvrvrvn fragile. Then, as the conceptual sketch shown in Fig. 6 for 

persists as the cut-on time-steps i max increases from f 0000 ° ' * ° 

nnnnn u „. , m, , , the potential barriers of the equilibrium states belonging 

to 20000 as shown m fug. 4. lhen, the activation energy * n 00 

a^i r t-i • 1 l ■ 1 r 1 ri / in\ \ to different classes, we may conclude that the states of 

AGi of Fq. (|9| is determined from the slope of log^r 1 - >) J 

i/n - -ua^i noi j a/~i as-* ir>- ft arc very easy to break up, and the states of Si are the 
vs. l/D, it yields AGq = 0.81 and AG\ = AG2 = 1.0, m- J J ^' 

, , , r,, /_ t=a 1 ht ,1 r . i most robust among the equilibrium states of three classes, 

dependent of the (7, GJ value. Moreover, the prefactor A; b H 

can be determined by the intersection between the straight As the state of collective behavior is of great inter- 
line of log (t^) vs. l/D and the vertical line of l/D = 0, ested, we then study its mean first-passage time for the 
we have Aq = 97.08, which is independent of the (7, G) system with the SF networks in a noisy environment. In 
value, and the other two prefactors A\ and A2 depend on the numerical calculations, we first assign a strongly disor- 
the network geometry. The 7-dependence may occur for der configuration to the system, then follow the stochastic 
Ag 1 when G is large, for example, we have As 1 = 10.11, LMR of Eq. © to generate a trajectory and to record the 
7.20, and 5.88 for G = 0.3 but 7 = 2.29, 2.55, and 2.82, time-steps for the first appearance of a state of So. Here, 
respectively; such dependence is not found for A$ 2 ■ On the the recorded time-steps is set as t max = 10 4 for the ab- 
other hand, both A\ and A2 are very sensitive to the G sence of the states of So at t — i max . We generate 10 6 
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trajectories totally for a given set of (7, C), and the mean nius equation. This allows us to quantify the robustness 
first-passage time, denoted as (to), is given as the av- of the equilibrium states through the values of the prefac- 
erage of the recorded time-steps over all trajectories. As tor and the activation energy in the Arrhcnius equation, 
the differences in (tq) caused by different 7 values are Our results indicate that the states of period-2 are very 
insignificant, we show the results as the plot of (tq) vs. fragile, and the trapped states are the most robust among 
D for 7 = 2.82 and C — 0.0066, 0.2016, and 0.3026, re- the three types. For systems in noisy environments, our 
spectively in Fig. 7. It is interesting to observe that the results obtained from the stochastic LMR indicate that 
(tq) value as a function of D is non-monotonic. As the there exists a range of noise for which, the mean first- 
D value increases, the (tq) value first decreases, reaches a passage time from strongly disorder states to the states 
minimum flat valley in the range of 0.07 < D < 0.22 with of system-wide coordination is the shortest. Thus, the ef- 
(tq) < 10, then increases abruptly at -D max — 0.22, and ficiency in reaching the state of collective behavior may 
the state of collective behavior becomes unreachable for be improved for systems with certain amount of noise. 
D > Dmax- Moreover, the role of clustering coefficient in As the distribution of subgroups in a social system can be 
the mean first-passage time is not very significant, as the characterized by the clustering coefficient, all these results 
results shown in Fig. 5 indicate, the difference in cluster- may provide wide applications in the study of collective 
ing coefficients is noticed by the mean first-passage time behaviors of social systems. 

for D < 0.10. Acknowledgement: We would like to thank Dr. Yutin 
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